A&A manuscript no. 

(will be inserted by hand later ] 



Your thesaurus codes are: 
06. (08.14.1; 02.04.1) 



ASTRONOMY 

AND 
ASTROPHYSICS 



Cooling neutron stars with localized protons 



o 
o 
o 

(N 

<: 



D.A. Baiko 1 and P. Haensel 2 

A.F. Ioffe Physical Technical Institute, Politekhnicheskaya 26, 194021, St. Petersburg, Russia 
e-mail: baiko@astro.ioffe.rssi.ru 
2 N. Copernicus Astronomical Center, Bartycka 18, 00-716, Warsaw, Poland 
e-mail: haensel@camk.edu.pl 

Accepted 20th January 2000 



> 

in : 
oo , 

o : 
o , 
o ■ 

^ : 

On; 

6 ' 
a: 



Abstract. We analyze cooling of neutron stars, assum- 
ing the presence of localized protons in the densest re- 
gion of their cores. Choosing a single threshold density 
for proton localization and adjusting neutron star mass, 
we reproduce the observational data on effective surface 
temperatures of Vela and PSR 0656+14, with or without 
an accreted hydrogen envelope. However, the presence of 
a tiny hydrogen envelope is mandatory, in this model, for 
reproducing the Gcminga data. 
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The model of dense neutron star matter with local- 
ized protons was proposed by Kutschera & Wojcik ( 1993Q . 
According to these authors, in high-density matter with 
low proton fraction the "zero-point" (Fermi) kinetic en- 
ergy of protons, estimated from the uncertainty principle, 
becomes small, compared to the energy of proton interac- 
tions with density waves of the neutron background. These 
density waves resemble hydrodynamic sound or long wave- 
length acoustic phonons in a solid. The coupling of pro- 
tons to the neutron density waves results in a significant 
increase of proton effective mass in close analogy with po- 
laron behaviour of a slow electron in a polar solid which 
also acquires a large effective mass because of the interac- 
tions with lattice phonons. 

Moreover, considering one proton in neutron matter, 
Kutschera & Wojcik have shown that above some critical 
density p\ p <; 4p (po = 2.8 x 10 14 g cm -3 ) the energy 
of a quantum state, in which a proton wave function is 
localized around some spatial point associated with the 
neutron density minimum, is lower than the energy of a 
state, where proton wave function is nonlocalized and the 
distribution of neutrons is uniform. Such a state also has 
an analog in solid state physics, the so-called small polaron 
state. It occurs when an electron interacts with a lattice 
so strongly that the latter deforms and this deformation 
traps the electron in a self-consistent manner. 

This result by Kutschera & Wojcik may be applied 
to neutron star cores, where one has finite admixture of 



protons (rather than one proton in neutron matter), if 
the proton fraction x p is so small that wave functions 
of neighbouring localized protons do not overlap. Small 
proton fractions in high den sity n eutron-star matter were 
obtained by Wiringa et al. ( 1988 ), who performed exten- 
sive many-body calculations of the ground state, based on 
variational method and assuming realistic nuclear hamil- 
tonian. In particular, their UV14+TNI model, which will 
be used further, predicts a proton fraction below 5% for 
densities p > 3po, and, eventually, complete disappearance 
of protons at about 7po- This result is in contradiction 
with some other calculations of equation of state (EOS) of 
neutron star matter [based on the relativistic mean-field 
approach, see Glendenning ( |1996| ), or performed within 
the Brueckner-Bethe-Goldstone many-body theory, e.g., 
Baldo et al. (1997)] which yield a proton fraction grow- 
ing monotonically with baryon density rib. The high den- 
sity behaviour of Xp(rib) is, therefore, subject to a signifi- 
cant theoretical uncertainty. Arguments in favor of x p (rib), 
which decreases and eventually vanishes at high nf,, were 
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The localization of protons (if confirmed) leads to dra- 
matic changes of basic kinetic properties of matter in the 
neutron star core. The reason is that the localized protons 
are very efficient scatterers of the main transport agents, 
electrons and neutrons. In effect, in the low-temperature 
regime, the thermal conductivity k, for instance, behaves 
as oc T, instead of the conventional T _1 dependence when 
scattering is provided by degenerate particles. As shown 
by Baiko & Haensel ( 1999| ) (hereafter Paper I), this ex- 
tends the time scale of thermal equilibration in the neu- 
tron star core by about two orders of magnitude indepen- 
dently of the state of neutrons (normal or superfluid). 

The neutrino emissivity of matter with localized pro- 
tons is also very non-standard. First of all, the con- 
ventional beta-processes which involve charged currents 
(direct or modified Urea processes) are forbidden since 
at given baryon density the minimum of energy corre- 
sponds to a specific proton fraction, while any process 
which changes this fraction would require rearrangement 
of the whole state and, consequently, much larger energy 
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than the thermal energy available, say, at T ^ 10 9 K. 
Thus, the only possible processes are those which involve 
neutral currents. These are bremsstrahlung of neutrino- 
antineutrino pairs in nn, np, and ep collisions. 

The kinetic properties of matter with localized protons 
were considered in Paper I. In particular, simple analytical 
formulae were derived for the electron and neutron ther- 
mal conductivities and shear viscosities, for the electron 
electrical conductivity as well as for the neutrino emis- 
sivity from np and ep bremsstrahlung. All these quanti- 
ties were calculated under two simplifying assumptions. 
Firstly, it was assumed that the localized protons did not 
exibit any correlations (i.e., behaved like impurities). Sec- 
ondly, the in-vacuum matrix elements of strong interac- 
tions were used. Both assumptions lead to an increase of 
the reaction rates by a factor of a few, since it is generally 
expected that the medium effects and correlation of scat- 
terers reduce the matrix elements of elementary scattering 
processes, at least for small momentum transfers. Thus the 
results of Paper I give the upper bounds to the neutrino 
emissivities and the lower bounds to the transport coeffi- 
cients. In other words, they represent the maximum effect 
that one can expect from the proton localization. 

The aim of the present paper is to apply the results of 
Paper I to modelling neutron star cooling in order to check 
if the hypothetical proton localization is compatible with 
observational data on the neutron star effective surface 
temperatures. To do that we have to specify the neutron 
star model. In the core (at densities p > p cx = 2.5 x 
10 14 g cm- 3 ), we take the UV14+TNI EOS by Wiringa 
et al. ( 1988| ), while in the crus t (be low p cx ) we take the 
EOS by Negele & Vautherin ( |1973| ) above neutron-drip 
density, and the EOS by Baym et al. ( ]1971 ) below neutron- 
drip density. The model is then determined by the only 
parameter, the central mass density p c , and is constructed 
by integrating the Tolman-Oppenheimer-Volkoff equation 
outward the center. 

We us e the cooling code described, e.g., in Levenfish et 
al. ( 1999] ) . The code includes explicitly effects of General 
Relativity (GR) and traces a steady thermal evolution of 
a spherically symmetric neutron star. Matter at densities 
P > Pb — 10 10 g cm -3 is assumed to be isothermal with 
constant temperature T?° (as measured at infinity). The 
local temperature, 7}, appears to be dependent on the ra- 
dial coordinate r due to the GR effect. The cooling is gov- 
erned by the thermal balance equation with energy losses 
due to neutrino emission from the core and photon emis- 
sion from the stellar surface. The local effective surface 
temperature T s is obtained from the local internal tem- 
perature, Tb, at p — pb using the formula by Potekhin et 
al. (1997). The latter formula fits the results of numerical 



simulations of heat transport through the crust. The for- 
mula depends also on the mass AM of the outer envelope 
made of light elements (light envelope) . Such an envelope 
may, in principle, reside on the neutron star surface and 
due to its higher thermal conductivity the envelope would 
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Fig. 1. Redshifted surface temperature versus neutron 
star age. Filled circles represent interpretation of observa- 
tions with the black-body model while open circles are the 
"atmospheric" interpretations. The threshold parameters 
for proton localization are x p q = 5% and p\ v = 4.3po- The 
solid curve corresponds to the stellar model with an en- 
velope of light elements, while the dashed and dot-dashed 
curves are calculated assuming no light envelope. The dot- 
ted curves are for the same star models (masses, presence 
or absence of the envelope) but without proton localiza- 
tion. 

decrease the thermal insulation of the inner stellar regions 
and modify the T s (Tb) relationship. The surface temper- 
ature T s °°, as measured by a distant observer, is related 
to T s as T s °° = T s y/1 - 2GM/Rc 2 , where M and R are, 
respectively, the mass and radius of the neutron star. 

For our purposes we have introduced a number of mod- 
ifications into the cooling code. All the changes are re- 
lated to matter at densities above p\ v . First of all, we 
have removed the proton contribution to the heat capacity. 
Then we have revised contributions to neutrino emissiv- 
ity. Above pi p , if neutrons are nonsuperfluid (which is the 
case in study), the main contribution to the neutrino emis- 
sivity is due to nn and np neutrino-pair bremsstrahlung. 
The rate of the nn process is not changed by the proton 
localization, while the emissivity of the np process reads 
(Paper I): 



<JK_ = 1-3 x 10" (i x ^ 



4n 



4/3 
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X 5* p Tg 



ergs s 



1 cm 3 , 



(1) 



where x p is the fraction of protons by number, n n is the 
neutron number density, m n and m* are the bare and 
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Fig. 2. Threshold density of localized protons versus mass 
of a neutron star for cooling model without envelope of 
light elements on the surface. The dotted and solid lines 
confine the parameter domains corresponding to the T s °° 
error bars of PSR 0656+14 and Vela, respectively. 

effective neutron masses, respectively; n — 0.16 fm~ 3 , 
and Tg is (local) temperature in units of 10 9 K. Finally, 
S\ p is a function of n n related to the total np cross-section 
in vacuum, 



S* p = 26.09- 



3.444- 



4nn 



(2) 



Typical cooling curves obtained within this model are 
presented in Fig. [l]. The effect of the proton localization 
can be easily seen by comparing the dotted curves, which 
represent the cooling without the proton localization, with 
curves of other types, which show a faster cooling due to 
the localization effect. The rate of cooling without the pro- 
ton localization is virtually independent of the mass: the 
dotted curves for 1.44 M Q and 1.48 M Q stars coincide. 
If the effect of the proton localization is taken into ac- 
count, the more massive star cools more rapidly because 
at fixed threshold density, pi p , the more efficient neutrino 
pair bremsstrahlung in neutron-localized proton collisions 
is operative in a larger region of the neutron star core 
leading to a larger overall neutrino luminosity. 

The results of our cooling simulations may be used 
for interpretation of thermal emission from the surface of 
cooling neutron stars of ages t ^ 10 4 yr, which is the typ- 
ical thermal equilibration time for neutron stars with lo- 
calized protons (Paper I). At present there are 4 sources, 
identified as isolated neutron stars, which show thermal 
soft X-ray emission and have characteristic ages t in ex- 
cess of 10 4 yrs. These are Vela, PSR 0656+14, Geminga 
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Table 1. The observational data on the 4 middle-aged 
isolated neutron stars. 



a > according to Lyne et al. (1996); 
[OFZ93] Ogelman et al. (|1993|) ; 

[PMC96] Possenti et al. (|l99^90% confidence level (c.L); 
[HW97] Halpern & Wang (|l997|) 90% c.L; 
[OF93] Ogelman & Fi nley ( |l993|) ; 

[PSZ96] Page et al. ( |1996|) , 90% c.L, for adopted lower 

limit of distance d = 300 pc; 

[ACPRT93] Anderson e t al. ( |1993D ; 

[MPM94] Meyer et al. ( |1994|) , typical value of the surface 
redshift factor, T s °°/T s = 0.8, were used. 



and PSR 1055-52. The data that we use are summarized 
in Table [|. The effective surface temperatures T s °° are de- 
termined from fitting the observed spectra of the sources 
in two different ways, either by black-body spectra (Tj°g 
in Table |l|) or by spectra obtained in realistic models of 
magnetized neutron star atmospheres (T^? in Table |l| for 
magnetized hydrogen atmosphere) . The inferred values of 
T s °° appear to be quite different. Radiation emerging from 
an atmosphere composed of light elements (hydrogen and 
helium) has harder spectrum than the black-body radi- 
ation with the same T s °° due to the decrease of opacity 
at higher energies (Romani 1987). Consequently, the ef- 
fective temperatures predicted by the atmosphere models 
are 2-3 times lower than the black-body ones which re- 
quire lower d/R (distance to radius) ratio. With the stel- 
lar radius fixed around the standard value of 10 km this 
translates into about 10 times smaller distances to the ob- 
jects. The atmospheres composed of iron, on the contrary, 
produce radiation with a softer spectrum which is more 
similar on average to that of a black body. 

The magnetic fields inferred from the spin-down rates 
are above 10 12 G for all 4 stars. Such fields must have 
significant effect on the atmospheric opacities and there- 
fore will modify the emerging spectrum (Shibanov et al. 
1992). In general, the magnetic field makes the spec- 
trum softer and somewhat elevates the effective temper- 
ature compared to the non-magnetized atmosphere case. 
However, no accurate magnetized atmosphere model have 
been developed so far for the relevant temperature range 
(T s < 10 6 K) where one should take into account eff ects o f 
atom motion on the opacity (e.g., Pavlov & Zavlin, 1998 ). 
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Fig. 3. Same as in Fig. || but for the model of a neutron 
star with envelope of light elements (of mass 10" 10 M) on 
the surface. The dashed line confines the (shaded) param- 
eter domain allowed for Geminga. 

These effects were neglected while obtaining the values 
reproduced in Table [j]. 

The surface temperatures T s °° are plotted in Fig. |l|. 
The filled and open circles correspond to black-body mod- 
els and simplified magnetized hydrogen atmosphere mod- 
els, respectively. As one concludes from Table || and the 
figure, the overall theoretical uncertainty in surface tem- 
peratures appears to be quite large; for instance T s °° ~ 
(2—6) x 10 5 K for Geminga. This does not allow one 
to draw any definite conclusion about the cooling sce- 
nario of the neutron star except that it probably rejects 
the "rapid" cooling (via direct Urea process unsuppressed 
by superfluidity in dense matter). We note that the sur- 
face temperatures obtained from the hydrogen atmosphere 
models are rather low and cannot be explained within the 
"standard" cooling scenario (the neutrino losses via mod- 
ified Urea process from non-superfluid matter). On the 
other hand, such a scenario can explain the black-body 
temperatures. Thus, the "atmospheric" temperatures (if 
the true temperatures turn out to be close to them) may 
provide more stringent test of the theory of neutron-star 
interior. For this reason let us adopt the "atmospheric" in- 
terpretation of observations and focus on lower error bars 
of Vela, PSR 0656+14, and Geminga. 

Having fixed EOS we are left with three parameters 
which influence the cooling. These are the central stellar 
density p c and the threshold values of mass density, pi p , 
and proton fraction, x p q. In addition, we can vary AM, 
the mass of the surface envelope made of light elements. 
Note, that the presence of the hydrogen atmosphere is 



not equivalent to the presence of a light envelope as the 
amount of hydrogen needed to modify the emerging spec- 
trum (1-10 g cm" 2 or ~ 10" 20 — 1O" 19 M ) is much lower 
than that required to change the T s (Tb) relationship. Let 
us fix the threshold proton fraction x p o = 5%, and con- 
sider two models: without envelope, Fig. 0, and with a 
light envelope of mass AM = 10~ 10 M, Fig. 0. The strips 
bounded by the lines of various types in FigsTg and ^| cor- 
respond to the domain of p\ v and M, for which cooling 
curves cross the error bar of a given source. Let us stress 
that Geminga (whose parameter domain is shaded in Fig. 
H) can be explained only assuming the light envelope. In 
the latter case and for realistic p\ v (presumably above 4po) 
the mass of Geminga should be above 1.5 M©. The other 
two sources can be explained either with or without the 
envelope but for fixed p\ p their masses in the model with 
envelope should be higher. Finally let us mention that 
the strips are not very sensitive to the envelope mass for 
AM > 10~ 10 M. 

Fig. [I] illustrates the cooling curves for the fixed 
threshold density p\ p — 4.3po- The solid curve represents 
the model with the light envelope and M = 1.61M© 
(p c = 5.4p )- It crosses the "atmospheric" error bars of 
Vela, PSR 0656+14 and Geminga simultaneously. The 
dashed and dot-dashed curves are envelope- free models. 
The masses should be lower, for instance, 1.44 and 1.48M© 
for Vela and PSR 0656+14, respectively. 

We have performed cooling simulations of neutron 
stars, with realistic equation of state, assuming localiza- 
tion of protons above some threshold density. These re- 
sults have been used for interpretation of effective surface 
temperatures of observed isolated, middle-aged neutron 
stars (t > 10 4 yr). We have shown that the available obser- 
vational data are consistent with the proton localization 
and can be reproduced if one chooses model parameters 
within the specific domains in the p\ p — M plane. Choosing 
a single threshold density for proton localization and ad- 
justing neutron star mass, we reproduce the observational 
data for Vela and PSR 0656+14, with or without an ac- 
creted, hydrogen envelope. However, the presence of the 
hydrogen envelope is required in this model for explaining 
the observations of Geminga. 

In conclusion, we note that the localization of protons 
is not necessary for the explanation of the available data 
on neutron star effective surface temperatures. Neverthe- 
less, as is seen from the Fig. [j], it changes significantly the 
cooling rate and, if its existence is confirmed, it represents 
an important effect to be taken into account in realistic 
neutron star cooling calculations. 
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